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I .  I ntroduct  f  on 

An  accurate  and  versatile  optical  model  of  stratospheric 
aerosols  has  important  applications  in  many  areas  of  geophysical 
research.  Such  models  are  indispensable  for  the  interpretation  of 
both  satellite  and  ground  based  remote  optical  sensors  that 
monitor  the  stratospheric  aerosol  layer.  They  can  also  provide  the 
basis  for  comparison  of  rather  diverse  types  of  measurements,  such 
as  evaluating  the  consistency  between  mass  and  scattering 
properties.  Stratospheric  aerosols  have  been  known  to  adversely 
affect  the  operation  of  remote  sensing  devices^  and  a  good  optical 
model  would  provide  the  basis  for  recognizing  and  possibly 
removing  the  degrading  effects.  In  addition  a  versatile  model 
should  be  able  to  provide  reliable  estimates  of  quantities  or 
parameters  that  have  not  been  measured  or  that  are  inconvenient  to 
measure.  Finally  it  might  be  observed  that  the  overall  success  of 
an  optical  model  is  a  direct  measure  of  the  intei — observational 
consistency  and  the  degree  to  which  the  aerosol  system  itself  is 
understood. 

Although  a  relatively  large  number  of  studies  have  been 
conducted  that  deal  with  stratospheric  aerosol  scattering 
properties,  only  a  few  efforts  have  been  directed  toward  a 
comprehensive  model  that  could  be  used  to  make  general  predictions 
and  form  the  basis  for  intei — instrumental  comparisons.  Shettle  and 
Fenn^  proposed  a  gamma  size  distribution  function  which  they 
considered  appropriate  for  optical  modeling.  The  variable 
parameters  In  this  function  allow  for  a  large  variety  of  size 
distribution  types.  However,  their  choice  of  parameters  was  not. 
necessarily  based  on  direct  size  distribution  measurements  in  the 
stratosphere. 

P innick  et  al^  proposed  a  simple  exponential  and  lognormal  size 
distribution  function  that  was  based  on  a  large  data  base  of 
directly  measured  submicrometer  size  aerosol.  The  model  was  used 
to  compare  particle  concentration  measurements  with  observed  1 idar 
backscatter ing.  It  was  also  used  to  estimate  the  dependence  of 
aerosol  extinction  and  backscatter  on  wavelength  and  to  calculate 
the  planetary  albedo  of  stratospheric  aerosols. 


Toon  and  Pollack^  proposed  a  global  stratospheric  aerosol  model 
based  on  a  zold  distribution  function.  Index  of  refraction  and 
total  optical  thickness  consistent  with  several  studies  available 
at  that  time.  (The  zold  distribution  has  been  shown  to  be 
equivalent  to  the  lognormal  dl strlbut lon^. )  This  model  was  used  In 
radiative  transfer  calculations  for  studies  In  long  term  climate 
effects . 

Russell  et  al^  made  an  effort  to  utilize  essentially  all  of  the 
previously  proposed  size  distributions  In  one  unified  approach. 

The  resulting  optical  model  provided  the  basis  for  comparing  a 
relatively  large  number  of  diverse  types  of  observations  made 
during  two  extensive  field  campaigns^’®. 

More  recently  Lenoble  and  Brogniez^  have  employed  bi modal  size 
distribution  functions  in  optical  model  calculations.  However, 
this  was  a  theoretical  study  and  the  size  distributions  themselves 
were  not  directly  matched  to  actual  measurements. 

All  of  the  above  optical  models  that  were  based  on  measured  size 
distributions  appear  to  be  relatively  successful.  However,  with 
the  possible  exception  of  the  work  of  Lenoble  and  Brognlez^  these 
models  were  applied  to  relatively  undisturbed  conditions  in  the 
stratosphere  and  it  was  not  possible  until  now  to  make  a  real  test 
of  their  general  applicability  to  a  wide  variety  of  natural 
conditions.  As  a  result  of  a  series  of  volcanic  eruptions 
beginning  with  Mt .  St.  Helens  in  May  1980  and  ending  with  El 
Chi  chon  in  April  1982,  the  stratospheric  aerosol  layer  has 
undergone  a  very  significant  change  in  both  concentration  and  size 
distribution*®.  The  present  period  (since  Jan.  1983)  is 
characterized  by  a  relatively  smooth  and  slow  decay  of  the  net 
disturbance.  This  stratospheric  perturbation,  which  may  not  occur 
again  within  the  next  100  years*®,  has  provided  an  unusual 
opportunity  to  further  test  and  advance  stratospheric  aerosol 
opt i ca 1  mode  1 s . 
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II.  Basic  Model  Development 


A.  Preliminary  Considerations 

The  particle  size  distr Ibut ion,shape  and  composition  are  usually 
at  the  heart  of  all  serious  optical  models.  Fortunately, 
stratospheric  aerosols  are  very  probably  essentially  spherical  and 
homogeneous  due  to  their  liquid  nature  (approximately  757.  solution 
of  sulfuric  acid)  and  therefore  MIe  scattering  theory  can  be 
applied  with  some  degree  of  confidence.  Recent  studies  of  the 
effect  of  composition  and  temperature  on  stratospheric  aerosol 
index  of  refraction  suggest  that  this  component  of  the  optical 
model  is  fairly  well  known  and  predictable^*.  Using  the 
calculation  procedure  of  Russell  and  Hamlll**  we  have  found  that 
for  all  of  our  soundings  (excluding  those  made  within  a  few  weeks 
after  a  volcanic  eruption)  the  Index  of  refraction  in  the  main 
part  of  the  aerosol  layer  (15  -  20  km)  varies  only  between  about 
1.44  and  1.45.  A  value  of  1.45  has  been  chosen  for  most  of  this 
work  since  it  Is  a  value  often  used  for  the  visible  wavelength 
region.  The  sensitivity  of  the  results  to  the  particular  choice  of 
index  of  refraction  will  be  discussed  In  a  later  section. 

The  primary  effort  in  developing  the  optical  model  here  is 
therefore  directed  toward  constructing  a  size  distribution  that  is 
consistent  with  a  number  of  well  determined  independent 
constraints.  The  authors'  past  experience  with  size  distribution 
measurements  Indicates  that  the  results  obtained  from  a  single 
method  or  investigator  should  never  be  accepted  Indiscriminately 
as  a  source  of  this  fundamentally  Imp>ortant  aerosol  feature. 
Optical  particle  counters,  for  example,  suffer  from  double  valued 
response  in  some  part  of  the  size  spectrum*^*  * ^ .  Impactor  type 
devices  are  frequently  used  to  sample  stratospheric  aerosols  and 
while  they  do  not  have  the  double  valued  response  problem,  some 
unknown  amount  of  evaporation  of  the  particles  may  take  place 
before  the  samples  are  processed.  In  addition,  there  is  some 
uncertainty  in  relating  the  distorted  shape  of  the  partic'les  on 
the  collecting  substrate  to  the  original  particle  diameter  In  the 
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free  stratosphere.  Size  dependent  collecting  efficiency  as  a 
funct ion  of  alt i tude  must  also  be  taken  i nto  account . 

In  this  study  we  generate  a  candidate  size  distribution 
consistent  with  our  extensive  data  set.  derived  from  a  balloon 
borne  system  of  optical  particle  counters.  Careful  consideration 
has  been  given  to  avoiding  particle  size  regions  where  the 
counters  experience  double  valued  response.  The  systematic 
accuracy  of  the  mass  and  scattering  properties  predicted  by  the 
candidate  distribution  is  then  tested  against  extensive  sets  of 
direct  measurements  made  by  other  independent  research  groups. 

This  procedure  differs  from  previous  studies  in  that  it  involves 
developing  a  consistent  size  distribution  from  the  constraints 
offered  by  several  extensive  data  sets  covering  a  wide  range  of 
stratospheric  aerosol  conditions.  In  this  approach  we  seek  to 
avoid  "tuning"  the  size  distribution  to  a  specific  stratospheric 
condition  or  to  the  results  of  a  few  specific  field  measurements. 

B.  Experimental  Procedure 

A  balloon  borne  system  of  three  individual  optical  particle 
counters  was  used  in  this  work  to  obtain  an  initial  estimate  of 
the  size  distribution  as  a  function  of  altitude.  One  counter  was 
dedicated  to  measuring  the  concentration  of  all  particles  greater 
than  about  0.01  micrometer  diameter,  which  for  the  stratosphere 
can  be  Interpreted  as  essentially  the  total  aerosol  concentration. 
Details  of  thi^  instrument  (  which  operates  as  a  high 
supersaturation  condensation  nuclei  counter)  have  been  presented 
e 1 sewhere^^ . 

The  second  optical  particle  counter  consists  of  our  standard  two 
channel  dustsonde^^,  which  determines  the  concentration  of 
particles  greater  than  0.15  and  0.25  micrometers  radius.  An 
estimate  of  the  slope  of  the  distribution  in  this  size  range  is 
given  by  the  ratio  of  the  counting  rates  In  channel  I  (radius  > 
0.15  wm  )  and  channel  II  (radius  >  0.25  um  )  and  will  hereafter  be 
referred  to  as  the  channel  count  ratio  I/I  I.  This  instrument  has 
been  In  use  for  more  than  20  years  and  is  known  to  produce  results 
consistent  with  a  wide  variety  of  other  types  of  Instrumentation 


for  background  aerosol  condl  t  fons^*®' The  past  success 
of  the  dustsonde  can  be  attributed  to  Its  relatively  accurate 
sizing  capabilities  (it  has  been  shown  that  this  counter  is  not 
affected  by  the  double  valued  response  region^^)  and  to  the  fact 
that  even  though  its  size  range  Is  very  limited,  it  Is  the  size 
range  most  relevant  to  visible  light  scattering  properties  for 
stratospheric  aerosol  background  conditions.  As  will  be 
demonstrated  below  however,  the  standard  dustsonde  does  not 
necessarily  cover  an  adequate  size  range  for  the  disturbed 
stratosphere. 

The  third  optical  particle  counter  used  in  the  compliment  of 
flight  instrumentation  is  a  dustsonde  modified  for  high  volume 
sampling  and  capable  of  measuring  particle  concentrations  greater 
than  four  different  radii:  0.25,  0.95,  1.20  and  1.80  urn.  The 
smallest  particle  channel  overlaps  with  the  standard  dustsonde 
measurement  and  we  require  agreement  between  the  two  instruments 
for  the  overall  sounding  to  be  considered  reliable.  The  light 
scattering  geometry  of  the  modified  dustsonde  is  essentially 
identical  to  that  of  the  standard  dustsonde  so  that  the  previous 
instrument  response  analysis  Is  for  the  most  part  still 
appi i cabi e^^ .  This  analysis  shows  that  the  1.20  and  1.80  ym 
channels  operate  well  out  of  the  double  valued  response  region 
while  the  0.95  ym  channel  is  very  close  to  it.  Considering  the 
fact  that  for  mono-dispersed  calibration  aerosols,  the  modified 
dustsonde  displays  a  somewhat  broader  pulse  height  spectrum  than 
the  standard  dustsonde,  a  revised  response  analysis  might  be  in 
order  and  could  put  the  0.95  ym  channel  in  the  double  valued 
response  region.  We  believe,  however,  that  the  1.2  and  1.8  ym 
channels  are  still  outside  of  the  double  valued  response  region. 
Thus  there  is  reason  to  exercise  a  good  deal  of  caution  in 
utilizing  the  0.95  ym  channel.  The  analysis  given  here  will 
proceed  along  two  paths:  one  in  which  the  Information  In  the 
potentially  suspect  channel  is  retained,  and  one  in  which  it  is 
ignored.  Certainly  it  is  not  Inappropriate  to  disregard  one 
channel  (only  a  loss  in  size  distribution  structure  would  result 
with  a  corresponding  smaller  probabllty  of  obtaining  a  useful  size 


df strf but  Ion)  but  rather  misleading  results  may  develop  by 
including  It.  From  the  beginning  it  was  recognized  that  there  may 
be  some  potential  problem  with  the  0.95  um  channel.  Our  reasoning 
in  choosing  the  location  of  the  four  size  channels  was  as  follows: 
the  0.25  um  size  provides  a  calibration  overlap  with  the  standard 
two  channel  dust sonde.  The  0.95  um  size  is  the  closest  value  of 
radius  to  0.25  um  that  still  has  a  chance  of  being  out  of  the 
double  valued  response  region.  The  1.20  um  size  was  considered  a 
backup  channel  to  the  0.95  um  channel  in  the  event  this  latter 
channel  displayed  suspicious  results.  The  fourth  channel  was 
chosen  for  as  large  a  particle  size  as  possible  consistent  with 
obtaining  a  workable  counting  rate  as  governed  by  the  sampling 
Flow  rate. 

C.  Data  Base 

The  six  channel  optical  particle?  counter  system  was  first  flown 
in  November  1980,  (after  the  disturbance  caused  by  the  Mt.  St. 
Helens  eruption)  but  regular  soundings  did  not  start  until  mid 
1982.  At  present  about  35  of  these  types  of  soundings  have  been 
conducted  from  Laramie.  This  paper  will  initially  focus  on  the 
data  obtained  after  January  1,  1983,  because  that  date  marks  the 
beginning  of  a  smooth  decay  mode  in  which  the  stratosphere  was 
relatively  well  mixed  but  still  of  a  considerably  different 
character  than  the  background  aerosol.  We  believe  that  this  is  an 
optimal  time  for  comparing  data  bases  generated  from  the  various 
observational  techniques  employed  at  mid-latitude  locations. 

D.  Size  Distributions 

From  the  above  considerations  it  is  seen  that  there  are  at  niost 
six  experimentally  determined  points  available  for  defining  the 
entire  integral  size  distribution  curve  and  that  a  rather  large 
gap  in  information  exists  between  0.25  and  about  1  um  radius.  The 
first  step  in  this  modeling  effort  centers  around  developing  a 
method  for  finding  a  "smooth"  curve  passing  through  the 
experimentally  determined  data  points.  By  smooth  we  mean  that  the 
first  and  second  derivatives  should  be  continuous.  This  guarantees 
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that  the  differential  size  distribution  as  well  as  its  slope  is 
continuous  -  a  seemingly  basic  requirement  for  a  realistic  size 
distribution.  In  addition,  it  will  be  required  that  this  smooth 
size  distribution  satisfy  the  constraints  imposed  by  other 
measurements  as  described  below.  The  reader  will  recognize  that 
this  procedure  could  lead  to  many  size  distributions  satisfying 
the  same  constraints  within  the  experimental  accuracy,  but  we 
proceed  on  the  working  assumption  and  expectation  that  all  such 
distributions  will  be  similar  and  will  produce  similar  results. 

Before  directly  proceeding  with  the  invention  of  new  size 
distributions  some  remarks  concerning  the  adequacy  of  previously 
used  size  distributions  might  be  in  order.  During  relatively 
undisturbed  periods  it  was  found  that  the  two  channel  standard 
dust sonde  data  could  be  used  to  define  a  simple  exponential  size 
distribution  that  would  adequately  account  for  the  scattering  and 
mass  properties  of  the  aerosol .  Some  refinement  was  made  to 
this  size  distribution  model  by  including  a  third  channel 
measurement  in  which  the  total  aerosol  concentration  was 
determined  with  a  condensation  nuclei  counter.  With  these  three 
channels  of  information  it  was  possible  to  define  a  unique 
lognormal  distribution  which  was  equally  satisfactory  in 
explaining  the  observed  aerosol  mass  and  scattering  properties  but 
treated  the  small  particle  region  more  realistically^*^'®.  For 
future  reference  these  early  models  were  used  to  calculate  the 
background  properties  where  illustrated  in  the  applicable  figures. 

The  exponential  and  lognormal  size  distributions  determined  from 
the  two  or  three  channel  data  after  the  1982  eruption  of  El 
Chi  chon  failed  quite  badly  in  that  they  were  not  consistent  with 
either  the  five  or  six  channel  size  data  and  were  not  consistent 
with  the  observed  aerosol  mass  and  scattering  properties.  Not 
surprisingly,  the  reason  for  the  failure  can  be  attributed  to  the 
fact  that  the  two  or  three  channel  size  data  cannot  simply  be 
extrapolated  to  correctly  account  for  the  concentration  of  the 
larger  particles.  Even  though  the  early  models  probably  did  not 
accurately  treat  the  large  particle  size  range,  they  were  still 
successful  because  the  concentrations  in  this  size  region  were  too 


small  to  effectively  contribute  to  the  aerosol  mass  and  scattering 
properties.  We  are  thus  forced  to  consider  some  new  functional 
forms  to  describe  the  perturbed  size  distribution  following  the  El 
Chi  chon  erupt  ion - 

An  exhaustive  discussion  dealing  with  the  merits  and 
deficiencies  of  the  various  size  distributions  tested  In  this 
study  would  be  quite  tedious  and  lengthy.  We  therefore  present 
only  a  summary  of  the  reasoning  leading  to  one  successful  choice. 

Some  success  has  been  achieved  by  fitting  the  six  channel  size 
data  to  a  two  mode  lognormal  size  distribution^^.  This  procedure 
seems  to  be  relatively  satisfactory  for  the  decay  period  following 
the  Mt -  St.  Helens  eruption,  (excluding  the  observations  of  the 
original  cloud  passing  over  Laramie)  but  not  for  the  period 
following  the  El  Chi  chon  eruption  because  it  generally 
overestimates  the  aerosol  mass  and  scattering  properties  and  gives 
a  suspiciously  narrow  mode  near  the  suspect  0.95  urn  radius 
channel.  At  first  this  result  seemed  plausible  because  it  was 
consistent  with  higher  resolution  size  distribution  measurements 
made  by  other  groups  during  the  same  time  period^®'^^.  At  present, 
however,  a  comparison  with  the  other  methods  Indicates  that  the 
concentration  implied  from  the  0.95  urn  channel  measurement  is 
considerably  (over  a  factor  of  10)  too  large.  The  explanation  for 
this  discrepancy  could  well  be  that  this  particle  channel  is 
influenced  by  the  double  valued  region  of  the  instrument  response 
as  discussed  above.  We  have  temporarily  abandoned  efforts  to  fit  a 
two  mode  lognormal  function  to  the  size  distribution  but  still 
recognize  that  the  possibility  exists  for  finding  a  viable  means 
of  utilizing  this  form. 

The  spline  curve  under  tension  is  a  useful,  but  brute  force, 
method  for  passing  a  smooth  curve  through  unequally  spaced  data 
points23.  By  adjusting  the  tension,  the  curve  changes  from  a  cubic 
spline  to  straight  line  segments  between  the  data  points.  The  fit 
can  be  made  in  either  linear  or  log  space  for  each  of  the  two 
variables  (concentration  and  radius  in  this  case)  which  gives  a 
total  choice  of  four  possibilities,  but  only  two  produce 
acceptable  results.  Curve  fits  to  quasi -power  law  distributions 


(as  might  be  encountered  fn  the  troposphere)  are  most  successful 
in  log- log  space.  On  the  other  hand,  quasi -lognormal  or 
exponential  curves  (as  past  experience  indicates  might  be  expected 
in  the  stratosphere)  are  best  fit  in  log  concentration  -  linear 
radius  space. 

The  size  distributions  obtained  with  spline  fits  under  tension 
utilizing  all  six  dust sonde  channels  generally  give  unsatisfactory 
results  by  overestimating  aerosol  mass  and  scattering  effects  and 
by  requiring  high  tensions  that  distort  the  distribution  curves. 
This  result,  along  with  the  above  concerns  expressed  for  the 
integrity  of  the  0.95  iim  channel,  has  led  us  to  abandon  its  use 
until  its  suspicious  behavior  can  be  satisfactorily  resolved. 

A  five  point  spline  fit  (neglecting  the  0.95  um  channel)  was 
also  tested.  In  this  case  satisfactory  results  could  be  obtained 
by  adjusting  the  value  of  the  tension.  However,  it  was  frequently 
necessary  to  use  relatively  high  values  of  tension  so  that  the 
differential  size  distribution  did  not  take  on  negative  values  in 
the  larger  size  region.  As  a  consequence  of  the  high  tensions  the 
size  distribution  became  unrealistically  distorted  in  the  smaller 
size  ranges.  This  behavior,  coupled  with  the  fact  that  there  was 
no  method  of  independently  selecting  the  correct  value  of  tension, 
led  us  to  abandon  the  approach  of  using  a  simple  spline  fit. 

One  of  the  most  successful  functional  fits  that  we  have  found  is 
a  combination  of  a  lognormal  and  spline  curve.  In  this  case  a 
lognormal  curve  is  uniquely  fit  to  the  first  three  size  channels 
(as  was  done  for  the  undisturbed  stratosphere)  and  truncated  at 
0.25  wm  radius.  A  spline  curve  under  tension  is  then  matched  to 
the  lognormal  curve  and  utilizes  the  remaining  data  points  except 
the  0.95  pm  channel.  The  resulting  size  distribution  has 
continuous  first  and  second  derivatives.  This  type  of  curve 
corrects  many  of  the  problems  associated  with  simple  spline  fit. 
The  lognormal  portion  of  the  curve  inherently  treats  the  small 
particle  range  of  the  distribution  in  a  more  realistic  fashion  and 
correctly  utilizes  the  total  measured  aerosol  count.  Furthermore, 
adjusting  the  tension  will  not  distort  or  affect  the  distribution 
in  the  small  particle  range  and  create  unrealistic  situations.  It 


should  also  be  noted  that  this  model  preserves  much  of  the 
functional  form  that  was  so  successful  fn  describing  the 
undisturbed  stratospheric  aerosol. 

In  evaluating  this  1 ognorma 1 -sp 1 i ne  (LNSP)  model,  calculations 
were  carried  out  in  both  log-log  and  log  concentration  -linear 
radius  space.  Similar  results  were  achieved  in  both  cases  when  a 
tension  of  approximately  2  was  used  for  the  log- log  space 
calculation.  The  results  of  the  log  concentrat i on- 1  inear  radius 
space  calculation  were  essentially  independent  of  tension  for 
values  between  zero  and  ten.  For  higher  tension  values  the 
distribution  curves  generally  became  unacceptably  distorted. 

Figure  1  illustrates  two  extreme  size  distributions  at  the  peak 
of  the  stratospheric  aerosol  layer  that  were  created  with  the  LNSP 
model  in  log  concentration  linear  radius  space  (  the  plot  is  in 
log  -  log  space  however).  As  can  be  seen  the  curves  are  not  very 
sensitive  to  the  value  of  tension  but  do  show  some  suspicious 
structure  for  the  highest  value  chosen.  Table  I  gives  the 
functional  form  and  parameters  describing  the  size  distributions 
shown  In  Figure  I  for  low  values  of  tension. 

The  results  of  the  calculations  for  the  LNSP  mode?  that  will  be 
illustrated  here  have  been  done  in  log  concentration  -  linear 
radius  space  because  they  are  essentially  independent  of  tension 
and  as  such  are  uniquely  determined  without  any  adjustable 
parameters.  This  would  not  be  the  case  if  the  calculation  were 
done  in  log  -  log  space  where  the  value  of  the  tension  can  have  an 
appreciable  affect  on  the  results,  and  an  independent  method  of 
selecting  the  appropriate  tension  would  need  to  be  available. 

E.  Size  Distribution  Constraints 

Using  the  LNSP  size  distribution  model  and  an  appropriate  value 
for  the  index  of  refraction,  the  total  aerosol  mass  concentration 
and  backscatter  cross  section  were  calculated  for  .5  km  altitude 
intervals  throughout  the  stratospheric  aerosol  layer.  These 
results  were  then  compared  to  the  direct  measurements  (at  the 
appropriate  altitude  )  as  reported  by  independent  observers.  Only 
those  independent  observations  were  used  for  which  we  were  able  to 
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obtain  a  significant  data  base  covering  the  smooth  and  weM  mixed 
decay  period  following  the  El  Chi  chon  eruption. 

Figure  2  shows  the  predicted  sulfate  mass  compared  to  the  actual 
measured  values.  The  calculations  were  made  for  a  particle 
specific  gravity  of  1.65  which  corresponds  to  approximately  a  757. 
sulfuric  acid  solution.  Each  value  was  computed  for  an  altitude 
corresponding  to  the  actual  measurements.  The  independent  mass 
values  were  determined  from  an  aircraft  borne  filter  sampling 
system  flown  near  Laramie  within  an  average  of  10  days  from  the 
corresponding  balloon  sound ing^^' 26.  xhe  data  covers  an 
altitude  range  of  13.7  to  19.8  km  and  a  time  period  extending  from 
April  1981  to  November  1983  when  the  high  altitude  mass  sampling 
program  was  unfortunately  suspended.  Figure  2  also  indicates  the 
degree  of  uncertainty  in  the  measurements  as  well  as  the  predicted 
values.  The  uncertainty  in  the  calculated  values  arise  from  the 
uncertainty  in  the  dust sonde  measurements  themselves,  which  has 
been  discussed  elsewhere^.  Considering  the  fact  that  there  are  no 
adjustable  parameters  Involved,  the  agreement  between  prediction 
and  measurement  is  quite  satisfactory  for  this  period  which  covers 
a  very  large  range  of  stratospheric  conditions  and  a  respectable 
altitude  interval .  We  have  therefore  Judged  the  LNSP  size 
distribution  model  as  being  acceptably  consistent  with  actual  mass 
measurements . 

Figure  3  shows  the  predicted  aerosol  backscatter  properties 
compared  to  available  lidar  measurements  made  at  0.6943  um 
wavelength  as  reported  by  the  NASA  Langley  Lidar  Group22.  por 
reasons  relating  to  availability  of  referenceabl e  data, 
comparisons  were  made  only  at  the  point  in  the  layer  where  the 
quantity  ’’backscatter  i  ng  ratio  -1"  was  a  maximum,  which  was  almost 
always  between  15  and  20  km.  The  comparison  covers  the  smooth 
decay  time  period  from  January  1983  to  January  1985  and  each  data 
point  corresponds  to  observations  made  during  the  same  month.  Also 
shown  in  Figure  3  i s  an  indication  of  the  uncertainties  in  the 
comparison  that  are  traceable  to  the  measurements  themselves.  The 
relatively  good  agreement  between  measurement  and  prediction 
indicates  that  the  LNSP  size  distribution  model  is  acceptably 


consistent  with  the  actual  aerosol  backscatter i ng  properties  at 
0.6943  urn  wavelength. 

Figure  4  illustrates  a  comparison  similar  to  that  shown  in 
Figure  3  but  corresponding  to  a  wavelength  of  1.06  um.  The  direct 
lidar  backscatter  measurements  were  made  in  Fukuoka  Japan  as 
reported  by  M.Hirono^^  and  cover  the  same  time  period  as  Figure  3. 
Even  though  a  small  systematic  difference  between  measurement  and 
prediction  may  be  indicated  in  Figure  4,  we  still  consider  that 
the  LNSP  size  distribution  model  is  acceptably  consistent  with 
direct  measurements  within  the  uncertainties  involved.  The  great 
distance  between  observation  sites  and  an  important  difference  in 
latitude  (Fukuoka : 33 . 6°N ,  Laramie:41 . 3°N)  makes  the  significance 
of  the  relatively  minor  discrepancies  rather  difficult  to  assess. 

In  our  procedure  for  developing  a  useful  size  distribution,  the 
LNSP  model  might  now  be  considered  acceptable  because  it  satisfies 
the  predetermined  constraints.  More  constraints  will  be  included 
in  this  process  if  and  when  additional  types  of  data  bases  become 
available.  As  previously  noted  the  LNSP  size  distribution  is  not 
unique  but  we  proceed  with  the  application  of  this  model  on  the 
working  tentative  assumption  that  it  will  not  produce  results 
significantly  different  from  another  model  conforming  to  the  same 
constraints.  More  work  in  the  future  may  be  needed  to  test  this 
assumpt i on. 

III.  Appi i cat 1 ons 

A.  Magnitude  of  Disturbance 

Before  continuing  with  strictly  optical  model  calculations,  the 
reader  should  be  aware  of  the  range  of  conditions  and  the  time 
period  for  which  this  model  applies.  The  peak  mass  mixing  ratio 
and  average  mass  radius  as  calculated  from  the  model  is  shown  as  a 
function  of  time  in  Figure  5.  The  figure  indicates  that  a  large 
perturbation  in  both  aerosol  mass  and  size  was  experienced  during 
the  study  period.  It  is  useful  to  note  that  only  after  the  El 
Chfchon  eruption  did  the  particle  size  change  s ign I f leant  1 y . 


Figure  5  also  shows  the  channel  ratio  I /I  I  and  fn  the  discussions 
below  we  will  show  that  this  parameter  can  be  indispensable  for 
the  interpretation  of  several  types  of  optical  measurements. 

B-  Interpretation  of  Lidar  Measurements 

It  has  been  suggested  that  the  Ildar  backscatter ing  to  mass 
ratio  (bks/mass)  could  be  relatively  constant  for  a  wide  variety 
of  size  distributions^®.  The  LNSP  model  provides  a  basis  for 
testing  this  potentially  simplifying  relationship  for  the  wide 
variety  of  size  distributions  encountered  during  the  study  period. 
Figure  6  shows  a  plot  of  bks/mass  as  a  function  of  I/II  at  the 
mixing  ratio  peak  for  a  wavelength  of  0.6943  um.  The  dashed  lines 
indicate  the  range  of  values  obtained  throughout  the  entire 
stratosphere  for  all  soundings.  It  is  clear  that  bks/mass  is  not 
strictly  constant  for  stratospheric  aerosols  as  hoped  but  the 
ratio  is  restricted  to  a  fairly  narrow  range.  It  is  also  clear 
that  a  much  better  estimate  of  bks/mass  could  be  made  If  the 
channel  ratio  I/II  were  available  from  dustsonde  measurements. 
Figure  6  provides  a  useful  basis  for  converting  lidar  backscatter 
data  to  mass  loading  for  all  altitudes  In  the  stratosphere. 

Figure  7  is  similar  to  Figure  6  but  the  calculation  was 
performed  for  a  wavelength  of  1.06  wm  and  an  index  of  refraction 
of  1.44  which  is  consistent  with  the  longer  wavelength.  In  this 
case  the  bks/mass  ratio  varies  within  wider  limits  and  it  is  more 
important  to  know  the  I/II  parameter  for  converting  the 
backscatter  measurement  to  a  mass  loading  estimate. 

The  model  also  provides  a  basis  for  comparing  the  consistency  of 
Ildar  measurements  made  at  different  wavelengths.  Figure  8  shows 
the  relation  between  the  calculated  backscatter  and  wavelength  for 
a  few  soundings  covering  the  range  of  aerosol  conditions.  The 
similarity  of  the  curves  suggests  that  it  would  be  relatively 
unproductive  to  attempt  to  derive  much  information  concerning  the 
size  distribution  of  stratospheric  aerosols  from  a  study  of  lidar 
backscatter  as  a  function  of  wavelength. 

B.  Optical  Depth 
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Using  an  approach  similar  to  that  of  the  previous  section,  the 
model  can  be  used  to  relate  ground  based  solar  extinction  or 
optical  depth  measurements  to  aerosol  mass  loading.  However  in 
this  case  the  measurements  usually  include  the  tropospheric 
component  and  as  such  cannot  be  directly  applied  to  assessing  the 
stratosphere  unless  the  tropospheric  optical  depth  is  known  to  be 
relatively  unimportant.  At  present  the  optical  model  is  useful  in 
estimating  the  contribution  of  the  stratospheric  optical  depth  to 
the  total  optical  depth.  This  type  of  study  will  be  continued  when 
large  data  sets  of  optical  depth  covering  the  relevant  time  period 
become  ava i 1 ab 1 e . 

Figure  9  illustrates  the  consistency  of  the  ratio  of  optical 
depth  (above  15  km)  to  the  column  mass  loading  (above  15  km) 
during  the  period  of  this  study.  It  would  appear  that 
stratospheric  optical  depths  could  be  converted  to  mass  loading 
values  fairly  accurately  for  column  integrated  channel  l/II  ratio 
values  greater  than  about  2.5.  For  lower  values  of  the  I /I  I  ratio, 
however,  some  knowledge  of  this  ratio  would  still  be  useful  in 
making  the  conversion. 

The  calculated  dependence  of  stratospheric  optical  depth  on 
wavelength  is  shown  in  Figure  10  for  several  representative 
periods  throughout  this  study.  Unlike  the  corresponding 
backscatter  calculations,  the  wavelength  dependence  of  optical 
depth  takes  on  a  completely  different  character  after  the  El 
Chi  chon  eruption.  Direct  observations^^  of  the  optical  depth  show 
very  similar  results  to  those  illustrated  in  Figure  10.  This 
change  in  the  relationship  between  the  optical  depth  and 
wavelength  might  be  considered  a  signature  of  the  El  Chi  chon 
eruption  in  the  stratospheric  aerosol  layer.  Furthermore,  these 
results  suggest  that  it  is  at  least  potentially  possible  to 
extract  some  useful  size  distribution  information  from 
measurements  of  atmospheric  extinction  of  solar  radiation  at  two 
or  more  wavelengths.  Such  studies  have  already  been  reported^®  but 
deal  with  size  distribution  types  not  necessarily  characteristic 
of  the  full  range  of  stratospheric  conditions. 
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C.  Interpretation  of  Satellite  Based  Solar  Extinction  Measurements 

Measurements  of  the  intensity  of  solar  radiation  traversing  the 
earth's  atmospheric  limb  from  a  satellite  during  spacecraft 
sunrise  and  sunset  represent  an  Important  and  useful  method  for 
remotely  sensing  and  monitoring  the  stratospheric  aerosol 
layer^»®.  An  accurate  optical  model  must  be  available  to  initially 
verify  the  operation  of  the  satellite  sensors  and  data  reduction 
algorithms,  and  later  to  relate  the  satellite  measurements  to 
other  parameters  such  as  particle  concentration,  mass,  and  1 idar 
backscatter.  In  this  application  the  LNSP  model  produces  results 
similar  to  those  shown  in  Figures  9  and  10,  but  specific 
applications  will  not  be  Illustrated  here  because  there  is 
unfortunately  no  mid- latitude  satellite  extinction  measurements 
available  covering  the  time  of  this  study. 

D.  Scattering  Phase  Function 

The  scattering  phase  function  Is  an  Important  aerosol  property 
In  many  types  of  radiation  transfer  calculations.  Figure  11 
illustrates  the  scattering  phase  function  predicted  from  the  LNSP 
model  at  the  peak  of  the  aerosol  layer  for  a  typical  sounding 
after  the  El  Chichon  eruption  and  is  very  similar  to  those 
reported  by  Volz^^.  Henyey  and  Greensteln^^  have  suggested  a 
useful  parameterization  of  these  types  of  phase  functions.  The 
bas 1 c  form  Is  g 1 ven  by 

5(0)  =  (l-g2)/4ff(l  +  g2  +  Zgcose  )3/2 

where  5(0)  Is  the  normalized  phase  function,  0  Is  the  scattering 
angle  and  g  is  an  adjustable  asymmetry  parameter.  This  function 
has  been  fit  to  the  model  prediction  shown  In  Figure  11  (dashed 
line)  and  as  can  be  seen,  the  agreement  Is  very  good  In  the 
forward  direction  but  relatively  poor  In  the  backward  direction. 
Nevertheless,  the  Henyey  -  Greenstein  function  Is  often  thought  of 
as  acceptably  representing  a  working  approximation  to  the  actual 


phase  function.  Figure  12  shows  the  variation  of  g  In  the  maximum 
of  the  aerosol  layer  during  the  study  period  as  calculated  from 
the  model.  As  shown  the  value  of  g  does  not  remain  strictly 
constant  but  a  relatively  accurate  value  could  be  estimated  from  a 
knowledge  of  the  channel  ratio  I/II. 

E.  LIdar  Backscatter  at  10.6  gm 

The  magnitude  of  aerosol  backscatter  near  10  um  wavelength  is  of 
critical  Importance  in  evaluating  the  feasibility  and  design  of 
the  WINDSAT  global  wind  monitoring  system^^.  Developing  a  global 
climatology  of  this  parameter  from  direct  measurements  would  be 
expensive  and  time  consuming  to  such  an  extent  that  it  may  well  be 
more  cost  effective  to  estimate  the  backscatter  from  other  types 
of  already  available  measurements  through  the  use  of  an  acceptable 
optical  model.  In  applying  the  LNSP  model  to  such  an  approach, 
comparisons  of  the  type  illustrated  In  Figures  2  and  3  should 
first  be  made  to  val Idate  the  model  for  this  wavelength  region  and 
test  the  consistency  of  the  observations  themselves.  Although  the 
data  exists  for  this  comparlson^^  the  study  has  not  yet  been 
successfully  completed,  but  ft  wll ?  eventually  offer  another 
important  constraint  for  atmospheric  optical  models. 

F.  Sensitivity  of  Scattering  properties  to  Index  of  Refraction 

The  LNSP  model  has  been  used  to  estimate  the  sensitivity  of 

aerosol  extinction  and  backscatter  to  variations  or  uncertainties 
in  index  of  refraction  for  the  range  of  size  distributions 
encountered  during  this  study  period.  It  was  found  that  for 
indexes  of  refraction  near  1.45  and  0.7  um  wavelength  the 
extinction  decreased  between  0  and  4%  for  a  change  of  -0.01  In  the 
index  of  refraction.  For  the  same  conditions  the  backscatter 
decreased  by  5  to  12X.  As  previously  discussed  the  range  of 
indexes  of  refraction  as  calculated  from  our  flight  data  using  the 
model  of  Russell  and  Hamill*^  is  about  1.44  to  about  1.45.  The 
extinction  and  optical  depth  measurements  are  therefore  not  very 
sensitive  to  inadvertent  fluctuations  In  the  Index  of  refraction. 


IV.  Conclusion 

It  would  appear  that  It  is  possible  to  develop  a  relatively 
satisfactory  optical  model  of  stratospheric  aerosols  from  a  rather 
low  resolution  size  distribution  measurement  supplemented  with  key 
constraints  derived  from  other  measurements.  The  specific  model 
employed  here  is  not  unique  and  further  work  should  be  done  to 
test  the  sensitivity  of  the  results  to  other  models  satisfying  the 
same  constraints-  A  more  satisfactory  functional  form  for  the  size 
distribution  should  be  found  so  that  the  model  itself  could  more 
easily  be  communicated  to  the  end  users.  At  present  our  entire 
data  base  must  be  accessed  to  make  calculations  covering  the  range 
of  stratospheric  conditions. 

This  work  was  primarily  supported  by  the  Air  Force  Geophysics 
Laboratory,  Hanscom  Air  Force  Base,  Bedford,  Mass.  Additional 
suppjort  was  received  from  NSF  and  NASA.  We  are  indebted  to  E.J. 
Mroz  for  allowing  us  to  use  his  unpublished  data. 
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Figure  1.  The  LNSP  model  size  distribution  for  two  extreme 
conditions  at  the  peak  of  the  stratospheric  aerosol  layer.  Each 
curve  was  constructed  for  tensions  of  .1.  3  and  30  as  marked.  Note 
that  there  Is  very  little  difference  between  tensions  .1  and  3, 
but  at  a  tension  value  of  30  the  curve  has  developed  more 
structure  than  might  be  considered  consistent  with  the  number  of 
data  points  defining  the  curve. 
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Figure  2.  A  comparison  of  the  LNSP  model  predicted  sulfate  aerosol 
mass  concentration  with  the  measured  values  derived  from  aircraft 
borne  filter  sampling.  The  magnitude  of  uncertainties  associated 
with  the  comparison  is  shown  for  one  point  only  but  is  typical  of 
all  points.  The  data  covers  the  time  period  from  April  1981  to 
November  1983  and  corresponds  to  discreet  altitudes  ranging  from 
13.7  to  19.8  km. 
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Figure  3.  A  comparison  of  the  LNSP  model  predicted  Ildar 
backscatter  at  0.6943  pm  wavelength  with  actual  measurements.  The 
magnitude  of  the  uncertainties  associated  with  the  comparison  Is 
shown  for  one  point  only  but  Is  typical  of  all  points.  The  data 
covers  the  time  period  from  January  1983  to  early  1985  and 
corresponds  to  the  peak  In  the  mixing  ratio  of  the  stratospheric 
aerosol  layer. 
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Figure  5.  The  time  development  of  the  peak  mixing  ratio  (ppbm)  and 
average  mass  radius  of  the  particles  In  the  stratospheric  aerosol 
layer  during  the  study  period  as  predicted  by  the  LNSP  model.  The 
horizontal  identifying  lines  at  the  left  Indicate  background 
values.  The  arrows  at  the  bottom  indicate  major  volcanic 
eruptions. 
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Figure  6.  The  predicted  dependence  of  the  backscatter  to  mass 
ratio  on  the  channel  ratio  according  to  the  LNSP  model.  The 
Individual  points  correspond  to  the  peak  of  the  stratospheric 
aerosol  mixing  ratio,  while  the  dash  lines  Indicate  the  range  of 
values  for  all  stratospheric  altitudes  up  to  about  30  km.  The 
uncertainties  In  each  point  due  to  uncertainties  relating  back  to 
the  measurement  Itself  Is  Insignificant  on  the  scale  of  the 
drawing. 
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Figure  7.  Same  as  Figure  6,  except  here  the  mode)  results  deal 
with  the  extinction  to  mass  ratio  at  1.06  um  wavelength. 
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Figure  8.  The  LNSP  model  predicted  wavelength  dependence  of  the 
backscatter  at  the  peak  of  the  stratospheric  aerosol  layer  for 
four  soundings  covering  the  range  of  conditions  observed  during 
the  study  period. 
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Figure  9.  Variability  of  the  optical  depth  /column  mass  ratio  as 
predicted  from  the  LNSP  model.  The  Integrated  channel  ratio  Is  the 
ratio  of  the  column  aerosol  concentration  (above  15  km)  In  channel 
I  (r>0.I5iim)  to  that  In  channel  II  (r>0.25wm). 
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10.  The  LNSP  model  predicted  wavelength  dependence  of 
depth  above  the  tropoMuse. 
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Figure  11.  A  typical  scattering  phase  function  at  the  peak  of  the 
aerosol  layer  after  the  El  Chi  chon  eruption  as  calculated  from  the 
LNSP  model.  The  dash  line  represents  the  Henyey-Greenstein  phase 
function  fit  to  the  model  results.  The  value  of  the  asymmetry 
parameter  g  Is  .73  In  this  case. 
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Figure  12.  The  LNSP  model  predicted  dependence  of  the  Henyey- 
Greensteln  asymmetry  factor  g  on  the  channel  ratio  I/I  I.  The 
Individual  points  were  derived  at  the  peak  of  the  stratospheric 
aerosol  layer*  and  the  dash  lines  Indicate  the  range  of  values 
obtained  throughout  the  stratosphere. 


